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1 Introduction 



. . .the ancients (as we are told by Pappus) esteemed the science of mechanics 
of greatest importance in the investigation of natural things, and the moderns, 
rejecting substantial forms and occult qualities, have endeavored to subject the 
phenomena of nature to the laws of mathematics . . . 
Isaac Newton, 1686 [30] 

Renormalized quantum electrodynamics is by far the most successful theory 
we have today. This very impressive fact, however, does not make the whole 
situation less strange. We start out from equations which do not make sense. 
We apply certain prescriptions to their solutions and end up with a power 
series of which we do not know that it makes sense. The first few terms of 
this series, however, give the best predictions we know. 
Res Jost, 1965 [16] 

In the more than 300 years that passed since Newton wrote this in his Principia Mathe- 
matical the moderns have been very successful at the endeavor to subject the phenomena 
of nature to the laws of mathematics - with exception of quantum field theory. As the 
second quote (which could have as well been written in 2002) shows, quantum field 
theory so far resisted a quantitative, mathematically rigorous foundation. 

In the present paper, an axiomatic approach is outlined that, I believe, provides foun- 
dations on which quantum field theory can be given a rigorous mathematical treatment. 
The present paper gives the elementary part and exhibits the connections to the tradi- 
tional settings. A deeper study of the consequences and use of the concepts presented 
here will be given elsewhere. 

In the new approach, each (classical or quantum) conservative physical system is char- 
acterized by two Hermitian quantities: a density and an action. A generalized Liouville 
equation defines the dynamics and implies Ehrenfest equations for expectations. 

For a classical (but not a quantum) field theory, the Ehrenfest equations in a symplectic 
Poisson algebra imply the traditional field equations by the stationary action principle. 
In particular, all traditional systems derivable from the stationary action principle can be 
modelled in our setting. In a similar way, one can get from suitable Lie-Poisson algebras 
relativistic and nonrelativistic Euler equations, Vlasov-Maxwell, Vlasov-Einstein, and 
Euler-Poincare equations. 

A new phase space quantization principle (generalizing the Wigner transform) allows 
the simple quantization of arbitrary Poisson algebras, with a good classical limit. 

A large class of Poincare invariant actions on spaces with a reducible representation of 
the Poincare group is exhibited. Since it is manifestly covariant but possesses a Hamil- 
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tonian nonrelativistic limit, it appears to be well-suited for phenomeno logical modeling 
of relativistic few-particle dynamics. 

In particular, we obtain a Lorentz-covariant phenomenological multiparticle quantum 
dynamics for electromagnetic and gravitational interaction which reduces in the non- 
relativistic limit to the traditional Hamiltonian multiparticle description with standard 
Newton and Coulomb forces. The key that allows us to overcome the traditional prob- 
lems in canonical quantization is the fact that we use the algebra of linear operators on 
a space of wave functions slightly bigger than traditional Fock spaces. 

For a quantum system, if the action is translation invariant, one can find pure states 
of given mass describing isolated systems in a rest frame by solving a constrained 
Schrodinger equation. This opens a constructive spectral approach to finding physical 
states both in relativistic quantum field theories and in phenomenological few-particle 
approximation. 

While I have already checked much of what is needed to get the many known results as 
consequences of the present setting, I am not yet completely sure about the adequacy 
of the new theory for all aspects of traditional field theory. Thus I'd like to apologize 
(as Newton did in the preface of [30]) and "heartily beg that what I have here done may 
be read with forbearance; and that my labors in a subject so difficult may be examined, 
not so much with the view to censure, as to remedy their defects. " 



2 Prelude: Covariant transmutation 

Do not imagine, any more than I can bring myself to imagine, that I should 
be right in undertaking so great and difficult a task. Remembering what 
I said at first about probability, I will do my best to give as probable an 
explanation as any other - or rather, more probable; and I will first go 
back to the beginning and try to speak of each thing and of all. 
Plato, ca. 367 B.C. [32] 

We begin with traditional nonrelativistic quantum mechanics of a multiparticle system 
and rewrite it in a formally covariant way that foreshadows the axiomatic setup devel- 
oped afterwards. (This section serves as a heuristic motivation only, without any claims 
to rigor.) 

Let if be a translation invariant and time-independent Hamiltonian, p the 3-momentum 
operator, the generator of the spatial translations, ip the energy eigenstate in the rest 
frame of a system with energy E > and total mass m. The Schrodinger equation gives 
Hip = Eip, and the condition that the system is in a rest frame says pip = 0. 
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To make these statements covariant, we extend wave functions by an additional argu- 
ment E. Then we may consider the energy E as an operator acting on functions ip(E, x) 
by multiplication with E, with time t = ihd/dE as conjugate operator. By introducing 
the operator 

L:=E-H, (1) 

we may write the Schrodinger equation together with the rest frame condition in the 
form 

Lip = 0, pip = 0. 

We now introduce a 4-momentum vector p = (^) , where p is related to the energy E 
by the relation 

E = p c-mc 2 . (2) 

On writing p 2 = p\ — p 2 (where p 2 = p • p) for the Lorentz square of a 4-vector, and 
applying a Lorentz transform, we see that this is equivalent with the condition 

Lip = 0, pip = kip 

for a pure momentum state ip with definite 4-momentum k > (in the forward cone, 
i.e., ko > |k|) and energy 

E = cVk? — mc 2 . 

A general pure state is a superposition ip = J dkipk of (unnormalized) momentum states 
with 4-momentum k, hence any nonzero ip with 

Lip = 0, i.e., ^GH:=kerL. 

A general mixed state is a mixture p = J dT[{a)ip a ip* ol of pure states ip a weighted by 

a nonnegative measure dn(a). Thus Lp = pL = 0. In particular, with the standard 
expectation 

(/>„:= trp/, 

of linear operators / on functions ip(E,x) we have 

[L,p]=0, (L) p = 0. (3) 

The equations (3) will be the starting point for our axiomatic setting. Slightly gener- 
alized, they will allow us to formulate not only nonrelativistic quantum mechanics, but 
also classical mechanics, classical field theory, and quantum field theory. 
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3 Axiomatic physics 



Das Streben nach Strenge zwingt uns eben zur Auffindung einfacherer Schlufi- 
weisen; auch bahnt es uns haufig den Weg zu Methoden, die entwickelungsfdhiger 
sind als die alten Methoden von geringerer Strenge. [...] 

Durch die Untersuchungen iiber die Grundlagen der Geometrie wird uns die 
Aufgabe nahe gelegt, nach diesem Vorbilde diejenigen Disziplinen axiomatisch zu 
behandeln, in denen schon heute die Mathematik eine hervorragende Rolle spielt; 
dies sind in erster Linie die Wahrscheinlichkeitsrechnung und die Mechanik. 
David Hilbert, 1900 [14] 

We now begin the axiomatic treatment; from now on, all concepts have a precise, un- 
ambiguous meaning. Here we concentrate on the conservative, (classical and quantum) 
mechanics part of Hilbert 's 6th problem, quoted above; for the probability part, viewed 
in the present context, see Neumaier [29]. In this paper, we only give the outlines 
and general flavor of the theory. A much more extensive version with full details, and a 
treatment of the dissipative case are in preparation. 

The quantities of interest are elements of a Euclidean Poisson algebra E containing 
the complex numbers as constants. Apart from an associative product (commutative 
only in the classical case) one has an involution * reducing on the constants to complex 
conjugation, a complex- valued integral J defined on a subalgebra IE of integrable 
quantities, and a Lie product (or bracket) ~i . The subalgebra -BE of bounded 
quantities consists of all / G E with /*/ < a 2 for some a G M. Quantities / with 
f* = f are called Hermitian. (For reasons given in Neumaier [29], we avoid using 
the customary word 'observables', and follow instead the International System of Units 
(SI) [43] in our terminology.) 

Apart from the standard rules for *-algebras and the linearity of the integral and the 
Lie product, one assumes the following axioms. (The product has priority over the Lie 
product, and both have priority over the integral. The partial order is defined by / > 
iff /* = / and J g*fg > for all g G IE, and the monotonic limit is defined by fi I 
iff, for every g G IE, the sequence (or net) f g* fig consists of real numbers converging 
monotonically to zero.) 

Axioms for a Euclidean *-algebra: 

(El) feBE,geIE g*,fg,gfem 

(E2) (.[//)• ./;/•• ffg = fgf 

(E3) fg*g>0 if g ± 

(E4) fg*fg = for all g G IE =^ / = (nondegeneracy) 

(E5) Jgf gi - Jf 9l - 0, Jgtf 9l ^0 

(E6) gi I =^ inf fgi = (Dini property) 
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Additional axioms for a Euclidean Poisson algebra: 

(pi) (f^ g y = r^ g * 

(P2) / n g — —g "i / (anticommutativity) 

(P3) / h (g n h) = (/ n gr) n ft, + # n (/ n h) ( Jacobi identity) 

(P4) / n gh— (f ~* g)h + g(f ~i ft) (Leibniz identity) 

(P5) /*/ = =^ / = (nondegeneracy) 

(P6) f ElE, g EE =>- f~\gem, 

(P7) ff~ig = if f£lE (partial integration) 

As a consequence, 

(P8) ff(gnh)=f(fng)h. 

Note that (E3) implies the Cauchy-Schwarz inequality 

KfgTifg) < Iff In n. 

which implies thast IE is contained in BE. 

The present definition of a Euclidean Poisson algebra is a modification of the concept of 
a Poisson algebra as discussed in Vaisman [45] and DA Silva & Weinstein [5] in that 
commutativity is dropped but integration requirements imposing a Euclidean structure 
are added. This modification enables us to treat classical and quantum physics on the 
same footing. (For a related attempt in this direction, see Landsman [21].) Moreover, 
we introduced the symbol ~i (an inverted stylized L, read 'Lie') to replace the Poisson 
bracket notation, which would be much more cumbersome if used extensively (as in as 
yet unpublished work). 

In the present, elementary paper, we make use only of some of the above axioms (mainly 
those not involving limits). However, as will be shown elsewhere, all are needed for the 
deeper analysis of our conceptual basis. 

To show that the axioms are rich in contents, we describe two basic realizations of them. 

The quantum Poisson algebra. Let H be a Euclidean (= pre-Hilbert) space. We 
define the commutator [/, g] := fg — gf, and let i := i/h with a positive real number h 
called Planck's constant. Then the algebra E = Lin H of continuous linear operators 
on HI is a Poisson algebra with quantum bracket 

f^g = t[f,g], 

and Euclidean with quantum integral 

J7 = W, 
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Integrable quantities are the operators / G E for which all gfh with g,/iGE are trace 
class. (This includes all operators of finite rank.) The axioms are easily verified. 

Nonrelativistic quantum mechanics. Nonrelativistic quantum physics is usually 
described by a rigged Hilbert space (see, e.g., Bohm [2]), if one wants to have direct 
access to the unbounded operators. Hence let H be a Euclidean space (the nuclear part 
of the rigged Hilbert space) with nuclear topology; we put EI = C°°(R, H ). For the 
standard position representation and po = ihd t /c, qo = ct, we have 

Pm -1 Qu = Vijv for/i, z/ = 0, 1,..., (4) 

with the metric 

77 = Diag(l,-l,...,-l). 
Restricted to E = Lin H , this gives the setting of traditional quantum mechanics. 

Nonrelativistic classical mechanics. As discussed, e.g., in Marsden & Ratiu [24], 
classical physics can be most conveniently described in terms of a Poisson manifold Q. 
Let {•, •} be the associated Poisson bracket on the algebra E := C°°(Q) of infinitely 
different iable complex-valued functions on Q. Then E = C°°(R 2 ,Eo) (defined as in 
Kriegl & Michor [19]) is a Poisson algebra with classical bracket 

= df dg df dg 
1 9 ' dtdE dEdt U,9i 

for / = f(t, E), g = g(t, E), and Euclidean with classical integral 

//= f dtdE f f(t,E) 
J Jn 

(where J Q is the Liouville measure). Integrable quantities are the Schwartz functions on 
E. (Thus integrability in the present sense is much stronger than Lebesgue integrability. 
This is due to our requirement (El) which implies that IE must be an ideal in E.) Again, 
the axioms are easily verified. With p — E/c, qo — ct and the standard symplectic 
Poisson bracket, we get again (4). 
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4 Physical systems 



. . . daft die iibliche Quantisierungsvorschrift sich durch 

eine andere Forderung ersetzen laftt [...J 

Die neue Auffassung ist verallgemeinerungsfdhig und 

ruhrt, wie ich glaube, sehr tief an das wahre Wesen der 

Quantenvorschriften. 

Erwin Schrodinger, 1926 [38] 

Motivated by the prelude, and consistent with the introductory remark in the seminal 
paper "Quantisierung als Eigenwertproblem" ("Quantization as eigenvalue problem") 
by Schrodinger [38], we generalize the Schrodinger picture of traditional quantum 
mechanics as follows. A physical system is characterized by a Hermitian density 
p G IE with p > 0. The density, or any set of parameters from which the density can 
be uniquely reconstructed by a well-defined recipe, is referred to as the state of the 
system. A physical system with density p defines expectations 

(/) == IPf = IfP- (5) 

The centralizer E(S) of a quantity S (or a vector of Lie commuting quantities) is the 
set of all quantities Lie commuting with (all components of) S, 

E(S) = {feE\S^f = 0}. 

Clearly, E(S) is again a Poisson algebra. For a quantity / G E(S), the conditional 
expectation at a fixed value s of S is defined by 

{f)s= s =(mS-s))/(5(S-s)), 

defined via a limit of integrable functions approaching the delta function. For example, 
if S is Hermitian with real spectrum then 

(/) s=0 = \im((S - ieY l f{S + ie)- 1 )/^ + e 2 )" 1 ). 

ej.0 

By construction, conditional expectations always satisfy (l)s=s — 1; they satisfy the 
axioms for an ensemble given in Neumaier [29]. 

Dynamical predictions are possible only in a system with a well-controlled environment. 
For a system in interaction with an arbitrary environment, the expectation satisfies a 
dynamics determined by the Ehrenfest equations 

(£>{/}) = for all / G BE (6) 
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with a forward derivation V, i.e, a continuous linear mapping V{-} : E — > E mapping 
bounded quantities to bounded quantities and satisfying 

v{fV = £>{/*}, v{f*f) > v{f*}f + f*v{f} 

for all / E E. Physical systems with the same forward derivation (but in general different 
densities) are said to follow the same physical law. Written in terms of the density, 
(6) becomes the generalized Liouville equation 

V*{p} = (7) 

with the Liouville operator V* defined (uniquely by (E4)) by 

fV*{p}f = JpV{f}. 

We shall discuss general physical systems and their (dissipative) properties elsewhere. 

Here we consider isolated systems only, where the physical law is characterized by a 
Hermitian action L G E which determines the forward derivation. A physical system 
with density p is called isolated (and p is called a conservative density) if 

(L) = JL P = 

and the generalized Liouville equation 

L n p = (8) 

holds. Since 

(L n /) = fp(L I /) = ftp I L)f = —J(L 1 p)f = 0, 
expectations in isolated systems satisfy the Ehrenfest equations 

(L n /) = for all / G IE. 

The axioms for a Euclidean expectation algebra imply that T>±{f} = ±L ~i / is a forward 
derivation for both signs; as discussed elsewhere, this reflects the reversible, conservative 
nature of isolated systems. 

If p is a conservative density and / G E(L) then 

P f := fpf* (9) 

is also a conservative density. Thus a large class of conservative densities can be con- 
structed from a single one if some quantities fi in the centralizer E(L) are known, since 
we may apply (9) with any polynomial constructed from the This generalizes the 
traditional construction of states from the vacuum by means of creation operators. It 
is applicable even where - such as for interacting quantum fields in 4 dimensions - no 
precise mathematical meaning can be given to the latter construction. 
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5 Hamiltonian systems 



Our axioms cover the traditional physics of Hamiltonian systems. The action corre- 
sponding to an arbitrary time-dependent Hamiltonian H(t) is defined as 

L = po-H(t), 

where p = E in the classical case and p = ihd t in the quantum case. In both cases, 

Po "> / = -/■ 

(Strictly speaking, the name 'action' fits tradition only for field theories. For multi- 
particle systems, the above expression for L is unrelated to traditional action principles. 
But applying the same machinery which gives the field equations of field theory to this 
unorthodox action happens to produce the correct multi-particle dynamics.) 

For conservative quantum systems, Lip = implies for L = p — H: 

P = -Po n P = ~{L + H)~t p— —H "i p, 

and we get the standard quantum Liouville equation 

ihp=[H,p] (10) 

for a conservative nonrelativistic quantum system with Hamiltonian H. The Ehrenfest 
equations reduce to their traditional form 

ih±(f) = (\f,n\), 

showing that expectations follow a deterministic law. For conservative classical systems, 
exactly the same derivation applies, and we get the classical Liouville equation 

p = {H,p}. (11) 



6 Pure states 

Pure states are the limiting situation (in a suitable completion of the space of integrable 
quantities) of densities extremal with respect to the natural order relation. They are 
of mathematical interest since any density can be written as a convex combination of 
pure states, and of physical interest for few-particle systems, where states can often be 
considered as approximately pure. (However, states at positive temperature are never 
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pure, and the decomposition into pure states is, in the quantum case, not unique. Thus 
pure states describe idealized situations only.) 

Pure classical states. Here extreme states are distributional limits of densities; the ex- 
pectations are algebra homomorphisms into C (i.e., characters of the algebra) satisfying 
(L) = 0. For nonrelativistic classical physics with phase space variable z, 

(f) = f(t,E,z), 

and the condition (L) = fixes the value of E to E — H(t, z). Hence we may assume / 
to be independent of E. 

Thus pure states of a classical nonrelativistic system are characterized by a pair (t, z) 
consisting of a time t and the phase space location z of the system at this time. The 
Ehrenfest equations reduce in the limit of pure states to the Hamiltonian dynamics 

f = {f,H}. 

Pure quantum states. Extreme states are limiting rank 1 densities 

p = $ G HP. 

The equation L ~i p = implies that ip is a generalized eigenvector of L. (See, e.g., 
Maurin [25] for a mathematical treatment in terms of nuclear spaces.) The condition 
(L) = then implies that the eigenvalue vanishes. Note that, since generalized eigen- 
vectors need not be in HI, not all expectations need to exist in a pure state; the latter 
are to be regarded only as idealized limits of physical states. 

Thus pure states of an isolated quantum system are characterized by a generalized 
Schrodinger equation 

Lip = o, ip e hi*. (12) 

We call solutions of (12) pure conservative quantum states. 

As discussed in the prelude, if the action L is translation invariant and p is the gen- 
erator of the translations, one can find pure conservative quantum states of definite 
4-momentum k by solving the equations. 

Poip = map, pip = kip. (13) 

In particular, pure states of mass m in a rest frame can be found by solving the eigenvalue 
problem 

p ip = map, pip = 0, Lip = 0. 
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This is a constrained Schrodinger equation, cf. Section 13 below. 



The pure conservative quantum states form a vector space EI cons on which the centralizer 
E(L) acts since / G E(L) and ^ e H cons imply L/ty = /L^ = /0 = 0. In the 
quantum case, / G E(L) iff / commutes with L; thus quantities in E(L) can be found, 
e.g., by solving the eigenvalue problem for L. Thus we can create from any particular 
conservative quantum state a large class of other conservative quantum states provided 
we know enough quantities commuting with L. 



We discuss here only boson fields. By using super Poisson algebras and super versions 
of all concepts, fermion fields can be handled in an analogous fashion. 

Let H := S^M 1 ' 3 ) be the algebra of infinitely differentiable, fast decaying Schwartz func- 
tions on Minkowski space M 1 ' 3 , and let V be a finite-dimensional symplectic space with 
symplectic form A. Then the field algebra E := C^iW® V*) of infinitely differentiable 
functions / of the field argument $ G H with d n f G C°°(M <g> V*, (H <g> \/)® n ) and at 
most polynomial growth is a Poisson algebra with 



With expolynomial functions (linear combinations of products of polynomials with the 
exponential of a negative definite, quadratic polynomial) as integrable functions, E is 
Euclidean with an integral definable via infinite-dimensional Gaussian measures. 

Pure classical field states. A pure state over the field algebra E = <g> V*) 

assigns to each / G E the value /(<E») at a particular field $ G H <g) V*. The Ehrenfest 
equations reduce in the limit of pure states over the field algebra to the equations 



7 Classical fields 




Li/ = for all / G E. 



Inserting the linear function / = a($), where 




into L "i / = we get 
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for a G H £g> V with compact support. Since A is nondegenerate, we conclude 

dL 



= for all x e M 1 ' 3 . 



d$(x) 

This is the traditional stationary action principle. In the current setting, it is not 
a postulate but a consequence of the Ehrenfest equations. (The equations for other 
choices of / are consequences of this.) 

To get the traditional field theories, we simply need to find the right symplectic structure 
for each type of field. The field components must appear in conjugate pairs, which we 
arrange to two conjugate vectors $ and $ c (in place of the single $ used before). Then 
adequate commutation relations are 

a($) n &($) = a($ c ) n 6($ c ) = 0, 
a($ c ) -i &($) = (a\b) := [ dxa(x) T b(x), 



where $ c = $* for complex fields (which come in complex conjugate pairs), while for 
real fields $ and $ c are independent. For real fields which have no conjugate partner 
in the Lagrangian, one adds additional conjugate partners to the algebra of quantities. 
These additional fields are - like gauge degrees of freedom - unobservable and do not 
affect the field equations for the original fields. 

Hence the present framework allows a consistent implementation of all classical field 
equations derivable from the stationary action principle. (Note: If we apply this to 
the electromagnetic 4-vector potential, we get, in contrast to the approach in canonical 
quantization, a conjugate 4-vector potential, with standard symplectic Lie bracket for 
each component!) 

By extending the above framework to Euclidean super Poisson algebras, one can also 
incorporate classical fermion fields. In particular, we can implement a classical version 
of the standard model, including gravitation within the present setting. 

If we use in place of symplectic Poisson algebras suitable Lie- Poison algebras, the Ehren- 
fest equations produce in the limit of pure classical states for appropriate actions both 
the relativistic [26] and nonrelativistic [27] Euler equations for perfect fluids and the 
Euler-Poincare equations [24]. Using suitable Lie- Poison algebras of functions of 
phase space fields, it is possible to define natural actions for which the Ehrenfest equa- 
tions produce in this way the Vlasov equations. In suitable tensor products one can 
then form actions that define Vlasov equations interacting with electromagnetic and/or 
gravitational fields, giving Vlasov-Maxwell equations (cf., e.g., [34]) and Vlasov- 
Einstein equations (cf., e.g., [1]). 

Details will be given elsewhere. 
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8 Phase space quantization 



There are many ways to quantize a classical system. From the point of view of being 
able to do analysis (i.e., error estimates), the mathematically most developed form 
is deformation quantization (see, e.g., Rieffel [37]), which deforms a commutative 
product into a Moyal product. In the following, we propose an alternative deformation 
approach which, instead, deforms the operators / e E by embedding E into Lin E, 
identifying / e E with the multiplication mapping g — > fg. This can be done with 
surprising ease. 

The superoperators Mj and Df defined by 

M f {g}:=fg, D f {g} := / n g 

belongs to Lin E. For / G E, we define the quantization / of / by 

- ih 

f := M f D f E Lin E. 

The expectations 

</> = </>-y (D f ) 

differ from those of / by a term of order 0(h), justifying an interpretation in terms of 
"deformation". In particular, we automatically have a good classical limit. 

To actually quantize a classical theory, one may choose a Lie algebra of relevant quanti- 
ties generating the Poisson algebra, quantizes its elements by the above rule, expresses 
the classical action as a suitably ordered polynomial expression in the generators, and 
uses as quantum action this expression with all generators replaced by their quantiza- 
tions. 

In general, the above recipe for phase space quantization gives an approximate Poisson 
isomorphism, up to 0(h) terms. But Lie subalgebras are mapped into (perhaps slightly 
bigger) Lie algebras, and one gets a true isomorphism for all embedded Heisenberg Lie 
algebras, i.e., Lie algebras where all Lie products are multiples of a central element 1. 

Quantization Theorem. If E is commutative then the quantum bracket 

A H B = l[A, B\ for A,B e Lin E 

satisfies, for /, g e E, 

^ ih 1 - - 

f^g = M f ^--Df^ g = -(Mf^ + f^g), 
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Any Lie subalgebra L of E defines a Lie algebra 



L = {M fng + h | f,g,heh} 

under the quantum bracket. If L is a Heisenberg Lie algebra then ~ : L — > L is a Lie 
isomorphism. 

The proof is not difficult but will be given elsewhere. 

In particular, for the standard symplectic Poisson algebra E = C°°(lR n x MI 1 ), phase 
space quantization amounts to using the reducible representation 

P = P- Y d n 1 = ( 1+ Y p 

of the canonical commutation rules on phase space functions instead of the traditional 
irreducible representation 

p = —ihd x , q = x 

on configuration space functions. It will be shown elsewhere that these representations 
are related by a Wigner transform (cf. Wigner [53]). 

By quantizing in phase space, one gives up irreducibility (and hence the description 
of a state by a unique density) but gains in simplicity. Perhaps this is comparable to 
the situation in gauge theory, where the description by gauge potentials introduces some 
arbitrariness with which one pays for the more elegant formulation of the field equations 
but which does not affect the observable consequences. 



9 Quantum field theory 

A good many physicists are now working on the problem 
of trying to set up a quantum field theory independently 
of any Hamiltonian. [. . .} 

I still think that in any future quantum theory there 
will have to be something corresponding to Hamiltonian 
theory, even if it is not in the same form as at present. 
Paul Dirac, 1964 [7] 

Actions for classical or quantum field theories are based on representations of a symmetry 
group and corresponding invariant actions. In any fundamental theory, the symmetry 
group must contain either the Galilei group (for nonrelativistic fields) or the Poincare 
group (for relativistic fields); if gravitation is involved, the symmetry group must also 
contain the group of all diffeomorphisms of some spacetime manifold. 
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Having a symmetry group is equivalent with having nonuniqueness in the description 
of a physical system. Different states providing equivalent descriptions (satisfying the 
same laws but in different coordinates) are commonly said to correspond to different 
choices of an inertial system. Changing the inertial system used to coordinatize a system 
changes the state and hence the expectations; for example, moving an inertial system 
O (illustrated by an observing intelligent robot) in time produces a change in observed 
expectations which O conceives of as the intrinsic dynamics of the environment, while 
moving (rotating or translating) the inertial system O in space produces a change in 
observed expectations which O conceives of as the illusion of the space moving around 
it caused by the motion of its moving head. We now formalize these considerations. 

Let L be the Lie algebra of the Galilei group, the Poincare group or any assumed 
symmetry group containing one of these groups, with Lie product ~i . Let p G L 1 ' 3 be 
the generator of the translation subgroup in the canonical basis. Let J be a Poisson 
representation of L in a Euclidean Poisson algebra E, defined by 

J (6) H J(5') = J{5 n 5') for all 5, 5' G L. (14) 

P := J(p) (taken componentwise) defines the (total) physical 4-momentum. A 
smooth change of the inertial system (modeling a virtual motion of the robots head) 
is described by an arbitrary continuously differentiable mapping 5 : [0, 1] — > L specify- 
ing the infinitesimal motions 5(t) £ L of the inertial systems at instant r G [0,1]. A 
corresponding assignment of densities p(r) G IE at instant r is called consistent if it 
satisfies the differential equation 

^p(r)=p(r)nj(5(r)). (15) 

In classical physics, this describes a canonical, in quantum physics a unitary transfor- 
mation representing a general element of the (connected part of the) symmetry group. 
In particular, an observer moving in space-time with uniform velocity u G R 1 ' 3 finds the 
density changing according to the covariant Liouville equation 

-^p(r)=p(r)nj( u -p). (16) 

Thus we have a covariant generalization of the nonrelativistic situation considered in 
Section 5. 

Since such a change of inertial systems should not affect the physics, we require that 
quantities (and in particular the action L, i.e., the physical law) are unaffected by these 
changes, and that an isolated system remains isolated. The former condition is simply 
the requirement that we base our setting on the Schrodinger picture, and the latter 
condition amounts to 

Lnp(r)=0 for all r (17) 
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whenever L ~i p(0) = 0. To analyze this condition, let p(0) be the density of an isolated 
system, and put 

e(r) := Lnp(r). 

Then e(0) = and 

^e(r) = ^(Lnp( r )) = Ln A p(r) = Ln(p(r)nJ(5(r))) 



so that 



= (Lnp(r))n J(5(r))+p(r)n(Ln J(5( r ))) 

^e(r) = e(r) n J(5(r)) + p(r) n (L n J(5(r))). (18) 

If (17) holds then e(r) vanishes identically, and this reduces to p(r) ~i (L ~i J(5(t))) = 0. 
The requirement that this holds for arbitrary densities and arbitrary smooth changes of 
the inertial system therefore demands that 

L n J (5) G CE for all 5 G L, (19) 

where CE denotes the Lie center of E, the algebra of quantities which Lie commute 
with all quantities. An action L satisfying (19) is called L-invariant. Conversely, if the 
action L is L-invariant, then (18) reduces to 

^(r) = e(T)i#)). 

Under conditions which guarantee the unique solvability of the initial value problem 
(18), we conclude that e(r) = for all r, proving (17). Thus the L-invariance of the 
action is essentially equivalent to the requirement that being isolated is a covariant 
concept. 

In particular, using in our setting a Poincare invariant action L defines a relativistic 
physical theory. As shown in Sections 4 and 6, we can use an arbitrary conservative 
density (resp. pure quantum state) and a set of quantities in the centralizer E(L) to 
construct a large class of conservative densities (resp. pure quantum states) as possible 
states of an isolated physical system with given action. 



Having phase space quantization as a universal generalization of the Wigner transform, 
we can use it to quantize the basic fields of any (Galilei or Poincare invariant) classical 
field theory. This gives well-defined mathematical definitions of the various (nonrela- 
tivistic or relativistic) quantum field theories in current physical usage. 

Using a Galilei invariant action one gets nonrelativistic field theory. As explained non- 
rigorously in many textbooks (e.g., Umezawa et al. [44]), nonrelativistic quantum field 
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theory is in principle equivalent to nonrelativistic quantum mechanics. Therefore, one 
uses for nonrelativistic problems field theory only to describe bulk matter, while scat- 
tering and bound state problems are handled with the Schrodinger equation. This is 
much simpler than solving the full operator dynamics of field theory. 

In relativistic quantum field theory, there has been in the past no analogue of the 
Schrodinger equation that could have been used for this purpose. Thus even simple 
scattering problems were formulated in a field theoretic language accessible to a per- 
turbative treatment, and bound state problems (see, e.g., Weinberg [47]) could be 
described only very indirectly through poles in the S-matrix. For the latter, there is 
no sound mathematical basis since in traditional quantum field theory, the S-matrix is 
only defined perturbatively in terms of a presumably divergent (Dyson [8]) asymptotic 
expansion, so that, mathematically, talking about its poles is nonsense. 

The results of the present paper show, however, that to each quantum (field or parti- 
cle) theory there is a corresponding constrained Schrodinger equation from which one 
can construct pure conservative quantum states with definite momentum in complete 
analogy to the nonrelativistic case, and without restriction to a particular symmetry 
group. (Mathematically, it is suspect if certain techniques work for a particular, highly 
nontrivial group but not for all groups. Already from this perspective one could see that 
something was missing from current quantum field theory!) 



10 Wight man axioms 

The quantum theory of fields never reached a stage where one could 
say with confidence that it was free from internal contradictions - nor 
the converse. In fact, the Main Problem [...J turned out to be [...] 
to show that the idealizations involved in the fundamental notions of 
the theory are incompatible in some physical sense, or to recast the 
theory in such a form that it provides a practical language for the 
description of elementary particle dynamics. 
R.F. Streater and A.S. Wightman, 1963 [42] 

Traditionally, mathematical physicists approach relativistic quantum field theory via 
an axiomatic approach discussed in detail by Streater & Wightman [42]. The 
Wightman axioms (Wightman [52]) are an interpretation of field theory not in terms 
of field equations but in terms of correlation functions. Relations to the Lagrangian 
approach have been lacking so far. But one would have such relations if one could 
combine tradition with the present formulation of quantum field theory. Thus one would 
like to realize the Wightman axioms by identifying the vacuum with a pure conservative 
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quantum state vpo with zero momentum, i.e., 

L^ = 0, P^o = 0, (20) 

and (in view of the remarks at the end of Section 6) Wightman field operators by suitable 
Hermitian quantities in the centralizer E(L). 

It is not clear whether the Wightman axioms describe correctly the structure of relativis- 
tic quantum states. Apart from generalized free fields, no realization of the Wightman 
axioms in 4-dimensional space-time is known (see, e.g., Rehren [33]), and there are 
no-go theorems - stating, for example, that there is no natural interaction picture [42, 
Theorem 4-16] - pointing to the possibility that these axioms are indeed too strong to 
describe realistic theories. 

To prove that the assumptions defining a Wightman field can (or cannot) be satisfied 
in the present context is therefore a highly nontrivial task. But at least it is embedded 
into a well-defined functional analytic context, where the Poincare representation is 
already fixed. This might make it tractable for systems like QED, which are close to 
nonrelativistic quantum mechanics. Therefore, one might be able to adapt the insights 
from nonrelativistic scattering theory (which provides a diagonalization of the action 
and hence full control over its centralizer) to the new situation. 

On the other hand, even without knowing the existence of Wightman fields (and even 
if one could prove that they do not exist), the setting presented here makes sense and 
defines for arbitrary actions a good quantum field theory, closely related to physical 
practice. In particular, one can try to generalize to the new constrained Schrodinger 
equations the supply of techniques available for ordinary Schrodinger equations, and in 
this way complement the current perturbative techniques of quantum field theory by 
techniques known from nonrelativistic quantum mechanics. A first step in this direction 
- the generalization of the projection formalism - has been done already; see Neumaier 
[28]. Work on scattering theory is under way. 

11 Phenomenological relativistic dynamics 

In spite of the acceptance of field theories as a matter of principle, 
most realistic dynamical calculations in nuclear physics, and many in 
particle physics, utilize the nonrelativistic Schrodinger equation. [...] 
Relativistic direct interaction theories of particles lie between local field 
theoretical models and nonrelativistic quantum mechanical models. 
B.D. Keistcr and W.N. Polyzou, 1991 [18] 

While fields are usually used to describe nature on a fundamental level, practical work 
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(especially for bound states and resonances) requires phenomenological few-particle 
equations, which are frequently related only loosely to underlying fields; see the ref- 
erences in the next section. It is therefore interesting to see that a variety of covariant 
phenomenological few-particle equations can be easily built in the present framework. 
We do this by using Poincare invariant actions on Hilbert spaces carrying a suitable 
Poincare representation without states of negative energy. 

The possible irreducible Poincare representations (modeling elementary particles) were 
classified by Wigner [54]. The representations of positive (relativistic) energy take 
their simplest form in momentum space; the momenta p are restricted to a mass shell 

£l(p) = {pe R 1 ' 3 | p 2 = p 2 , p > 0}, (21) 

the orbit of a 4-vector p under the Poincare group. It is possible to combine these 
irreducible Poincare representations in many ways to obtain reducible momentum space 
representations for few-particle systems. Traditionally (see, e.g., Weinberg [48] for 
the canonical field quantization approach and the review in Keister & Polyzou [18] 
for the direct relativistic Hamiltonian few-body approach), this is done by breaking 
the manifest invariance to a maximal subgroup of the Poincare group, with all the 
awkwardness this entails. 

The key that allows us to preserve a manifestly covariant formalism, thus overcoming 
the traditional problems in canonical quantization, is the fact that we use as algebra 
of quantities the linear operators on a space of wave functions slightly bigger than 
traditional Fock spaces. This is done in the following by adding a velocity vector u as a 
dynamical parameter, which allows us to deform the bare mass shell p 2 = (mc) 2 (where 
c is the speed of light) to p 2 = (mu) 2 , which in turn permits the conservation of total 
4-momentum in interactions. 

A phenomenological realization of a system of N massive scalar particles with rest 
masses m 1 , . . . , m N > and charges Q 1 ,. . . , Q N is now realized by wave functions 

ij) = ip(u,p 1:N ) = ^(u,p\ ...,p N ) 

whose coordinates are a global 4-velocity vector u with < u <E IR 1,3 and the particle 
4-momentum vectors p a in the dynamic mass shells Vl(m a u) whose scale depends 
on u. The total 4-momentum ^p a is required to be parallel to the 4-velocity u. Thus 
the space of wave functions is 

H = C°°(n N ), (22) 

where Q N is the set of all tuples 

(u,p 1 -- N ) = (u,p\...,p N ) 
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with 

p a E tt(m a u) for a — 1, ... ,N, 0<u \\ ^p a . 

The (not everywhere defined) Hermitian inner product - from which a Hilbert space can 
be constructed by completing the space of vectors of finite norm - is given by 

0*V := J dm duDp 1 . ..Dp N d(mu - ^p a )0(u,p 1:JV )^(u,p 1:JV ), 



where 



Dp = dp 5(p* - {muf) = ^= dp (23) 

2po 2^(mw) 2 + p 2 

is the invariant measure on a dynamic mass shell Q(mu). The one-particle operators 

are defined as 

J(f) Q ,m ,p ,M a ), 

a 

where the diagonal operator / = f(u,Q,m,p,M) is a function of 4- velocity u, charge 
Q, mass m, 4-momentum p and 4-angular momentum 

d 

M:=pA— (24) 
op 

with components 

=P^ Pu 



O fV rv 1 

OPu OPu 

and the superscript a indicates application to the coordinates of the ath particle. (Note 
that the global 4- velocity u carries no superscript; it is shared by all particles.) Since the 
p M are the Poincare generators of translation in the direction of the //-axis and the M^ u 
are the standard generators of the Lorentz transformations, it is easy to see that the total 
4-momentum J(p) and total 4-angular momentum J(M) define a representation of 
the Poincare group without negative energy states. In the terminology of DlRAC [6], it is 
a representation in the point form. (It shares this property with the representations of 
RuiJGROK [36] which are based on Lippmann-Schwinger equations. But his translation 
generators are much more complicated than the present ones.) 

On the space (22), one can now define actions of the form 

L = L - V, (25) 

where the kinetic action L is a Poincare invariant one-particle operator, and the 
interaction V is a Poincare invariant integral operator. 
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12 Poincare invariant multiparticle interactions 



For scalar particles, the simplest covariant kinetic action is 

w(? ^ M>2 _ c2)) , (26) 

with a constant c > 0, the speed of light. However, more complicated covariant 
formulas with rational or analytic dependence on m and p 2 are admissible, too, if they 
vanish for p 2 = (mc) 2 and nowhere else. In this case, the generalized Schrodinger 
equation Lip = implies for noninteracting particles, where V = 0, the relation u 2 = c 2 , 
forcing the dynamic mass shells to equal the bare mass shells. 

To construct a versatile class of Poincare invariant interactions, we first note that the 
vector 

— p ■ u + J Op ■ u) 2 — p 2 u 2 + (mu 2 ) 2 

p m :=p + u w - (27) 

u z 

is in the dynamic mass shell Q(mu). Indeed, it suffices by covariance to check the case 
where u = 0; then uq > 0, u 2 = u^, p ■ u = poUo, 

iPm)o =Po + u — 2 = VP + {muo) 2 > 0, 

u 

and since p m = p, we find p 2 m = (m-u ) 2 = (mu) 2 . Thus the mapping p — > p m (the 
dependence on u is not written explicitly) is a nonlinear projection to the dynamic mass 
shell VlOmu). 

The simplest choice for a nontrivial interaction is a sum of pair interactions, 

V = J2 ya \ (28) 

a<b 

where V ab = V ba acts on the coordinates of particles a and b as 

(V ab iP)(u,p a ,p b ) = Jdq 5(u ■ q)U ab (qMu, (p a + q) m *, (p b - q) m >), (29) 

where the projections are to be taken with respect to the common 4-velocity argument 
u, and U ab (q) is also allowed to depend on mass, momentum and charge of the particles a 
and b. The delta function removes a redundancy in the projections, which do no change 
if a multiple of u is added to q. The construction is such that V ab is automatically 
translation invariant. In particular, if all U ab are Hermitian and Lorentz invariant then 
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V and hence the action (25) is Hermitian and Poincare invariant. For example, this is 
the case in pair potentials of the form 



jjab 



(q) = Re 



(3 s {m a m b ) 2 + (I3 v m a m b + aQ a Q b )p a ■ p b + (3 T {p a • p b f 
m a m b (q 2 + ie) 



(30) 



where the limit e { is to be taken in (29) to regularize the potential near q = 0. 
These potentials describe relativistic electromagnetic and gravitational forces; the cou- 
pling constants a and Ps,Pv,Pt determine the strength of the electromagnetic and 
the scalar, vector, and tensor gravitational interaction, respectively. (This will be 
justified in the next section by considering the nonrelativistic limit.) By making these 
coupling constants g-dependent (running coupling constants), one can also account 
covariantly for phenomenological self-energy contributions; cf. the discussion in Peskin 
& Schroeder [31, pp. 252-255]. 

Note that after Fourier transform into spacetime, we get - in contrast to field theories 
- a nonlocal (but still Poincare invariant) action. 

This basic setting can be extended in various ways. Particles with positive spin or with 
internal symmetries are easily accommodated, especially when using the representations 
discussed in Weinberg [48, 49]. (They are of course equivalent to Wigner's represen- 
tations but computationally more tractable.) Particles with positive integral spin are 
handled in exactly the same way, except that the wave functions have additional indices, 
the angular momentum gets an additional intrinsic spin term operating on these indices, 
and the inner product has a slightly different form. It is easy to specify L-invariant in- 
teraction terms similar to (30) for particles with positive spin and for particles with 
inner symmetries (and corresponding matrix- valued charges Q a ); but such interactions 
are now also restricted by Clebsch-Gordan rules (cf. Weinberg [50]). 

Fermion particles with half-integral spin are handled similarly, using spinor components 
in the wave functions and kinetic actions such as 



generalize the Dirac equations to the multiparticle case. Details about the handling of 
spin will be given elsewhere. 

Massless particles are handled in the same way, except that the kinetic part of the action 
is absent, since these particles never go off-shell in our phenomenological setting. 



L Q = J(p ■ 7c — mc 2 ). 



(31) 



The resulting constrained Schrodinger equations 



J(p ■ 7c — mc 2 )ip = Vip, pip = kip 
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For indistinguishable particles, symmetrization and antisymmetrization can be done in 
the standard way. Different kinds of particles are handled by adding to the sum of their 
self-actions another interaction. Few-particle systems in which the particle number is 
not conserved can be modelled by using a direct sum of Hilbert spaces of the type (22) 
and covariant interactions changing the particle number. For example, we may model 
the emission and absorption of a photon of momentum p by a massive scalar particle of 
charge Q a with the Hermitian and Poincare invariant interaction proportional to 

(V^)(u,p,p a ) = F(p)^^(u, (p a + p) ma ), 

(V^)(u,p a ) = [ dp5(p 2 )F(p)^-^(u,p,(p a -p) m «), 
J m a 

where the form factor F(p) is an arbitrary covariant scalar C°°-function formed from 
p, p a and u. (Note that the photon wave function has additional vector components, 
with respect to which the inner product • is taken.) Previous covariant few-particle 
models could not handle this situation (Keister & Polyzou [18, p. 392]). 

In a multiparticle system, one can model in the same way the interactions correspond- 
ing to Feynman diagrams with a single vertex of degree 3, and in a similar way also 
interactions corresponding to more complex vertices. Note that because momentum 
is conserved and all particle energies are positive, particles cannot be created from a 
vacuum state (with particles), nor can particles be annihilated without creating (or 
preserving) at least one particle. Thus the phenomenological approach does not have 
the problems which field theories have with the presence of an interacting ('fluctuating') 
vacuum. 

We see that the possibilities for the new action-based relativistic models fully match 
(and even exceed) the freedom available for nonrelativistic Hamiltonian systems. Since 
they are manifestly Poincare invariant, they are much simpler than various relativis- 
tic Hamiltonian models that have been constructed in the past (see, e.g., the review 
in Keister & Polyzou [18] for nuclear physics, Crater et al. [4] for QED, and 
RuiJGROK [36] for a Lippmann-Schwinger based model), but have the same advantages 
as the latter: consistency with relativity theory, tractable few-body calculations, easy 
treatment of bound states, resonances, and particle production, and easy fit to paramet- 
ric models. In addition, they can be used to give phenomenological models of quantum 
systems in which the particle number is not preserved, or the spin is > 1. 

In time, such action-based relativistic models may therefore replace the many nonrel- 
ativistic (e.g., Isgur [15], Karl [17]), semirelativistic (e.g., Lucha et al. [22, 23]) 
and relativistic (e.g., Keister & Polyzou [18]) Hamiltonian approximations, and ap- 
proximations based on Bethe-Salpeter equations (e.g., Kummer & Modrich [20]) or 
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Dyson-Schwinger equations (e.g., Roberts & Williams [35]) now in vogue for the 
phenomenological description of quarks, mesons, baryons, and other relativistic matter. 
Since our phenomenological actions are easily made manifestly symmetric under the full 
symmetry group of a system, it may also give more workable low energy effective the- 
ories for the standard model, such as chiral perturbation theory (e.g., ECKER [10, 11]) 
or quantum hadrodynamics (e.g., Serot [39], Serot & WALECKA [39, 40]). 

The relation between the above action-based relativistic multiparticle models and the 
field-theoretic models discussed earlier is not clear at present. It is expected that the 
projection techniques from Neumaier [28] relate the field theories from Section 9 to 
corresponding effective iV-particle theories modeled as in the present section. On the 
other hand, it is also conceivable that the field theories should rather be regarded as 
limits of iV-particle theories in the thermodynamic limit iV — > oo. There are indications 
that this might be the case for QED (since radiation phenomena are always dissipative) 
and for gravitation (since black hole thermodynamics changes pure states to mixed 
states, cf. Wald [46, pp. 180-185]; the traditional coupling to a hydrodynamic model 
is also meaningful only in a thermodynamic limit). 



13 Constrained Schrodinger equations 

States of fixed total 4-momentum J(p) can be obtained by solving (13). With a Lorentz 
boost, we may transform the system to a rest frame; the resulting constraint J(p) = 
can be imposed kinematically by restricting the 4- velocity to u = 0. Since c and J(mc 2 ) 
are constants, the wave function is an eigenstate of the rest frame energy J(p c — mc 2 ) 
(a shifted relativistic energy p^c, introduced in analogy to the prelude), and we are left 
with the (still rotation invariant) constrained Schrodinger equations 

ifj = 6(u)i/> , Lxjj = 0, J(p c - mc 2 )ilj = Eip, (32) 

the relativistic analogue of the nonrelativistic multiparticle Schrodinger equation af- 
ter separation of the motion of the center of mass. Thus our phenomenological ap- 
proach is a covariant version of the situation in the prelude: The mass shells form 3- 
dimensional manifolds, and the momenta p a can be considered as relativistic analogues 
of 3-momentum vectors. Since u = 0, the 4-velocity contributes only one additional 
degree of freedom uq, which replaces the energy degree of freedom of the nonrelativistic 
situation. Thus, in contrast to the realizations of quantum field theory discussed above, 
to traditional Bethe-Salpeter equations, and to proper time based relativistic multipar- 
ticle dynamics (see, e.g., Fanchi [12]), there are no superfluous degrees of freedom, but 
the treatment is still manifestly covariant. 
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The delta function in the interaction (29) forces q = 0. Dropping the redundant 
coordinates u = 0, p = a/ (mU ) 2 + p 2 and q = from the notation, the interaction 
can be written as the 3- dimensional integral 

(V ab ij)(u ,p a ,p b ) = c- 1 J dqU ab (q)^(u ,p a + q,p b -q); (33) 

the prefactor comes from the delta function in (29). If we now Fourier transform in 
space to get the position representation, 

^Kx a ,x 6 ) = J rfp a rfpV pa - xa e tp6 - x V( Mo ,p a ,p 6 ), 

we find 

V°^p{u , x a , x 6 ) = U^(x b - x a )^K, x a , x 6 ) 
with the spatial potential 

U^(r) = c' 1 J dq e tqr f/ ab (q). (34) 

This looks like a nonrelativistic formula, but the covariant nature of the model is visible 
in the form (29) of U ab (q) and also shows in the constraint nature of (32). Compared to 
the nonrelativistic case, this is now a general linear eigenvalue problem for the eigenvalue 
E, and its solution is slightly more demanding. But numerical methods are available; 
see, e.g., Golub & van Loan [13]. 



The nonrelativistic limit. To deepen the analogy, we give a rough, heuristic deriva- 
tion of the nonrelativistic limit c — > oo; it would be interesting to have a rigorous 
version of this from which one can obtain error bounds. The equation Lip = can 
be written (for bosons) as J(|m(w 2 — c 2 ))ip = Vtp. For small potential energies, 
V J(m)c 2 and small spatial momenta, p 2 (mc) 2 , this gives u 2 = c 2 + 0(1), 
hence p 2 = (mu) 2 = (mc) 2 + 0(1) and p = a/ (mc) 2 + p 2 = mc + 0(c~ l ). Therefore, 

Pl ~ 2 ( ™ C)2 = (po - mc)^p = (po - mc)(c + 0(c' 1 )) = p c - mc 2 + 0(c- 2 ), 

L = J 



V 2m J V 2m J \2m) 



w J(p c - mc 2 ) - J(p 2 /2m) = E - J(p 2 /2m) 

with the rest frame energy E = J(poC — mc 2 ). Thus, the constraint Schrodinger equation 
reduces in the nonrelativistic limit to the standard Schrodinger equation for a multipar- 
ticle system with Hamiltonian 

H = J(p 2 /2m) + ^f^(x 6 - x a ). 

a<b 
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Arbitrary local pair interactions can be obtained in the nonrelativistic limit by choosing 
U ab appropriately (and in a non- unique way). For larger kinetic energies, the potential 
in position space acquires additional, nonlocal terms (that can be approximated using 
derivatives in the interaction). Thus we have a flexible covariant theory with a good 
nonrelativistic limit. 

In particular, from the covariant potential (28) with pair interactions of the form (30), 
we recover in the nonrelativistic limit the standard nonrelativistic multiparticle 
dynamics in the presence of electromagnetic and gravitational forces. 



14 Golden opportunities 

Behind it all is surely an idea so simple, so beautiful, that when 
- in a decade, a century, or a millennium - we grasp it, we 
will all say to each other, how could it have been otherwise? 
John Archibald Wheeler, 1987 [51] 

Eine mathematische Theorie ist nicht eher als vollkommen 
anzusehen, als bis du sie so Mar gemacht hast, daft du sie dem 
ersten Manne erkldren konntest, den du auf der Strafte triffst. 
David Hilbcrt, 1900 [14] 

I do not know whether the perfection requested by Hilbert can be achieved in deep 
theories. But, having discovered the unexpected beauty of the present approach, I hope 
that the insights presented will contribute to the perfection of quantum field theory. 

In 1972, Freeman Dyson [9] gave a lecture called "Missed opportunities", where he 
talked " about the contribution that mathematics ought to have made" to physics "but 
did nof . I believe the present contribution widely opens the door for mathematicians to 
contribute to quantum field theory, and creates golden opportunities for those interested 
in mathematical physics. 

The present setting gives a mathematically consistent point of view from which to study 
the laws of physics, which complements the point of view taken by past history. On 
the new basis, it is likely that scientists will resolve in the near future the most ba- 
sic challenges current theoretical physics poses to mathematicians and mathematical 
physicists: 

• the existence of QED and derivation of its properties, 

• bound states and resonances in quantum field theories, 
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• a unified quantum field theory of all forces of nature, 

• the existence and mass gap in quantum Yang-Mills theory - one of seven Clay 
millenium prize problems [3], a golden opportunity in the most concrete sense. 
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the results presented here. 

I also want to thank the maintainers of (and the contributors to) the Los Alamos Na- 
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read more thoroughly while searching for the path to success. 

I'd like to thank Dr. Hermann Schichl (Wien) for many discussions on various pieces of 
the puzzle that helped me to clarify my thoughts. Thanks also to Prof. Peter Michor 
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to Prof. Walter Thirring (Wien) for his treatise on mathematical physics which I used 
over and over again, to Prof. Gerhard Ecker (Wien) for useful discussions on quantum 
field theory long ago, and to Prof. Hartmann Romer (Freiburg), who introduced me 
many years ago to the idea that an elementary particle 'is' [41, p. 149] an irreducible 
unitary representation of the Poincare group. 
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